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ABSTRACT 

Some gamma-ray burst (GRB) afterglow light curves show significant variability, which 
often includes episodes of rebrightening. Such temporal variability had been attributed in sev- 
eral cases to large fluctuations in the external density, or density "bumps". Here we carefully 
examine the effect of a sharp increase in the external density on the afterglow light curve by 
considering, for the first time, a full treatment of both the hydrodynamic evolution and the ra- 
diation in this scenario. To this end we develop a semi-analytic model for the light curve and 
carry out several elaborate numerical simulations using a one dimensional hydrodynamic code 
together with a synchrotron radiation code. Two spherically symmetric cases are explored in 
detail - a density jump in a uniform external medium, and a wind termination shock. The ef- 
fect of density clumps is also constrained. Contrary to previous works, we find that even a very 
sharp (modeled as a step function) and large (by a factor of a 1) increase in the external 
density does not produce sharp features in the light curve, and cannot account for signifi- 
cant temporal variability in GRB afterglows. For a wind termination shock, the light curve 
smoothly transitions between the asymptotic power laws over about one decade in time, and 
there is no rebrightening in the optical or X-rays that could serve as a clear observational sig- 
nature. For a sharp jump in a uniform density profile we find that the maximal deviation Aa^ax 
of the temporal decay index a from its asymptotic value (at early and late times), is bounded 
(e.g, Ackmax < 0.4 for a = 10); Actmax slowly increases with a, converging to Actmax ~ 1 at 
very large a values. Therefore, no optical rebrightening is expected in this case as well. In 
the X-rays, while the asymptotic flux is unaffected by the density jump, the fluctuations in a 
Sire found to be comparable to those in the optical. Finally, we discuss the implications of our 
results for the origin of the observed fluctuations in several GRB afterglows. 
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1 INTRODUCTION 

Gamma-ray bursts (GRBs) are produced by a relativistic out- 
flow from a compact source. The outflow sweeps up the external 
medium and drives a strong relativistic shock into it. Eventually 
the outflow is decelerated (by pdV work across the contact discon- 
tinuity that separates the ejecta and the shocked external medium) 
and most of the kinetic energy is transferred to the shocked external 
medium (for recent reviews see Piran 2005; Meszaros 2006). The 
shocked external medium produces long lived afterglow emission 
that is detected in the X-rays, optical, and radio for days, weeks, 
and months, respectively, after the GRB. The afterglow emission 
is thought to be predominantly synchrotron radiation. This is sup- 



ported both by the broad band spectrum and by the detection of 
linear polarization at the level of a few percen t in the optica l (or 
near infrared) afterglow of several GRBs (see ICovindl2003l. and 
references therein). Inverse-Compton scattering of the synchrotron 
photo ns might dominate the observed flux in the X-rays in some 
cases jPanaitescu & Kumail200CtlSari & Esinl200l|:lHarrison et alJ 
I2OOII) . 

In the pie-Swift era the best monitoring of GRB afterglow 
light curves was, by far, in the optical. Most afterglow light 
curves showed a smooth power law decay iStanek et al,. 199^ 
iLaursen. & Stanekll2003l iGorosabel et"ai]E00 g), and often also a 
smooth achromatic transition to a steeper power law decay that 
is attributed to the outflow being collimated into a narrow jet 
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Rhoadslll999l. ISari Piran & Halpemlfl999h . It has been argued 
Wang & Loebll20(XI) that the smoothness of the afterglow light 
curve is directly related to (and thus enables to constrain) the 
smoothness of the external density field. Nevertheless, some op- 
tical afterglows have shown significant temporal variability, with 
strong deviations from the more typical sm ooth power law be- 
havior. The bes t examples are GRB s 021004 j Pandev et ai]|2002t 
Foxetal. 2003; Bersieretal. 2003) and 030329 ( Matheson et al. 
2003; Sato et al. 2003; Uemura et al. 2004; Liokin et al. 2004). 

Possible explanations for such temporal variability in GRB 
afterglow light curves include variations in the external density 
iWang&Loeb 2000; Lazzati et al. 2002; Nakar, Piran & Granot 
|2003D, or in the energy of the afterglow shock. The lat- 
ter includes energy injection by "refreshed shocks" - slower 
shells of ejecta that catch up with the a fterglow shock 
on long time scales i Rees & Meszaros '1998^, 'Kumar & Pirar/ 
boOOa; Sari & Meszaros 2000; Ramirez-Ruiz, Merloni & Rees 
I2OOII: iGranot. Naka r & Piran 200^) or a "patchy shell" - angu- 
lar inhomogeneities within t he outflow jK umar & Piran] |2(^(^(^^]|: 
[ Nakar. Piran & Granoll Eoosl iHevl & PemiJl2003t iNakar & OreJ 
l2004h . Another possible cause for variability in the afterglow light 
curve, although it is expected to be quite rare, is microlensing by an 
intervening star in a galaxy that happens to be close to our line of 
sight. GRB 00030 IC exhibited an achromatic bump in its optical 
to NIR light curve s that peaked after ~ 4 days iSagar et'al]l2OO0l : 
iBerger et al."200o') which had been interpreted as such a microlens- 
ing event ^Gar navich. Loeb & Stanek 2000; Gaudi, Grano t & Loebl 
l200lUBaltz & Hui 2005), although other interpretati ons such as a 
bump in the external density have also been suggested iBerger et alJ 

Recent observations by Swift have found flares in the early X- 
ray afterglows of many GRBs tBurrow s et al. 2005; Nousek et alJ 
l200d iFalcone et aLll2006l; lO'brien et alJl200d) whi ch are proba- 
bly due to late time activity of the central source ^Nousek et alj 
[2OO6; Zhang et al. 2006). Early optical variability also appears to 
be more common than previously thought (e.g.. .Stanek et all20o9) . 
although it is not yet clear if it is caused by similar mechanisms as 
the late time optical variability that had been detected before Swift. 

The most natural forms of variations in the external density 
are either clumps on top of a smooth background density distri- 
bution, or a global abrupt change in density with radius. The latter 
can be, e.g., the termination sho ck of the wind from the massive star 
progenitor of a long-soft GRB JWiiersll200lh . Such a stellar wind 
environment may have a richer structure and can include an abrupt 
increase in density with radius at the contact discontinuity between 
shocked wind from two different evolutionary stages of the pro- 
genitor star, as well as clump s that are f ormed due to Rayleigh- 
Taylor instability (Ramirez-Ru iz et alJ200 5 ; Eldridge et al. 2006 ) ' . 
Density clumps with a mild density contrast may also be formed 
due to turbulence in the external medium. Furthermore, the ex- 
ternal density p rofile is expected to vary between different pro- 
genitor models iFrver. Rockefeller & Youn3l2006h . The variabil- 



' iRamirez-Ruiz et all l2005t) find that the clump formation may also in- 
volve the Vishniac instability, and once such clumps are formed they stand 
a reasonable chance to survive until the time of the core collapse of the 
progenitor star. 



ity that had been observed in optical afterg lows was attributed 
to density clumps in the ex ternal medium iLazzati et alj |200'3 ; 
iNakar. Piran & Granot 2003). The expected observational signa- 
tures of the afterglow shock running into the wind termination 
shock of the massive star progenitor have also been considered 
<Wiiersl200lHRamirez-RuizetalJ200lll2005l ; lEldridgeetalJ200d 
IPe'er & Wiiersll2006h . In all these cases it has been argued that 
there would be a clear observational signature in the foiTn of a 
rebrightening in the afterglow light curve, before approaching the 
new shallower decay slope corresponding to the uniform density of 
the shocked wind. 



Here we revisit the effect of density fluctuations on the af- 
terglow light curve by solving in detail the case of a spherically 
symmetric external density with a single density jump (by a fac- 
tor of a > 1) at some radius Ro, while the density at smaller and 
larger radii is a (generally different) power law in radius. This is 
done by constructing a semi-analytic model for the observed flux 
due to synchrotron emission at different power law segments of 
the spectrum and by carrying out numerical simulations. The semi- 
analytic model takes into account the effect of the reverse shock on 
the hydrodynamics and on the emissivity, as well as the effect of the 
spherical geometry on the arrival time of photons to the observer. 
Being semi-analytic, however, this model uses some approxima- 
tions for the hydrodynamic evolution and the resulting radiation. 
Therefore, we also perform numerical simulations in which the 
light curves are calculated using a hydrodynamic-l-radiation numer- 
ical code. This code self-consistently calculates the radiation and 
evolves the electron distribution in every fluid element, which cor- 
responds to a computational cell of the one dimensional Lagrangian 
hydrodynamic code. The results of this code are used to obtain light 
curves in cases of special interest and near spectral break frequen- 
cies, as well as to verify the quality of the semi-analytic model, 
which is found to agree well with the numerical results. These cal- 
culations are much more accurate than those presented in previous 
works, and our results are significantly different. In all cases we find 
a very smooth transition to the new asymptotic power law, with no 
rebrightening for an initially decaying light curve. 

In ^we develop the semi-analytic model. First, in 32.11 a 
simple analytic model is constructed for the hydrodynamics, which 
agrees very well with our numerical results. Then, in 32.2l we con- 
struct a semi-analytic model for the observed flux density. Specific 
case studies (^ are then analyzed in detail, for a wind termination 
shock ( 33. H and for a spherical density jump in a uniform medium 
( 33. 2> . The light curves for these cases are also calculated using 
the numerical code (described in Appendix IaI. The effect of prox- 
imity to a break frequency around the time of the density jump is 
investigated in 33.31 and our main conclusions are found to remain 
valid also in the vicinity of the break frequencies, ^is devoted 
to a discussion of the expected observational signatures of density 
clumps on top of a smooth underlying external density distribution. 
Such density clumps are found to have a very weak observational 
signature which would be very hard to detect. In ^we discuss the 
implications of our results for the origin of the observed fluctua- 
tions in several GRB afterglows. Our conclusions are discussed in 
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2 SEMI-ANALYTIC MODEL FOR A SPHERICAL JUMP 
IN THE EXTERNAL DENSITY 

In this section we model a spherical relativistic blast wave that 
propagates into a power-law external density profile (p 

ext — A/" 

with A: < 3) which has a single sharp density jump (by a factor 
a> 1) at some radius r = Ro. The power law index, k, of the exter- 
nal density is allowed to be different at r < {ko) and al r > Ro 
(Ml 

The hydrodynamic evolution, as well as the resulting contribu- 
tion to the light curves, can be roughly separated into three phases 
corresponding to the following ranges of the radius R of the forward 
shock: (i) aiR < Rq the blast-wave follows a self-similar evolution 
( Blandford & McKee 1976, BM hereafter), (ii) atR = Ro a reverse 
shock forms which crosses most of the shell of previously shocked 
material (that had been swept up at r < Rq) alR = Ri, while the for- 
ward shock continues ahead of the density jump but with a reduced 
Lorentz factor, (iii) alR> R[ the forward shock relaxes into a new 
self-similar evolution corresponding to the new density profile at 
r > Ro- In the following we first approximate the hydrodynamic 
evolution of the different shocks during these three phases and then 
calculate the resulting light curves. 



2.1 Hydrodynamics 

Consider a spherical ultra-relativistic blast wave (identified with 
the afterglow shock), which is well described by the self-similar 
BM solution atR< Rq, that propagates into the following external 
density profile: 



Pext '• 



'AQr"" r<Ro 



Air-'' r>RQ 



(1) 



The amplitude of the density jump, i.e. the factor by which the den- 
sity increases at r = Rq, is given by 



lim 



Pext[(l + e)/^o] Al k k 

= —Rn 



-opext[(l-e)/?o] Ao 



(2) 



and is assumed to be larger than unity. The afterglow shock en- 
counters the jump in the external density profile at a lab frame time 
ro = [1 + 1 /2(4-ko)rl]Ro /c ^ Ro /c, where c is the speed of light ,r4 
is the Lorentz factor of the shock front just before it encounters the 
density bump at r = Ro, and the corresponding Lorentz factor of the 
fluid just behind the shock is denoted by 74 = r4 /y/l. At f < ?o there 
are three regions: the region behind the afterglow shock (subscript 
'4') is described by the BM solution with (A,k) = (Ao.ko), and the 
two regions of cold unperturbed external medium (subscripts '0' 
and '1' air <Ro and r > Rq, respectively). 

When the afterglow shock encounters the jump in the exter- 
nal density a reverse shock is driven into the hot BM shell, while 
a forward shock propagates into the cold higher density external 
medium at r > i^o- At this stage region '0' no longer exists, but 
two new regions are formed so that altogether there are four re- 
gions: (i) the cold unperturbed external medium ahead of the for- 
ward shock with a density pext =Air"*', (ii) the shocked external 



medium originating from r > Ro, (iii) the portion of the BM shell 
that has been shocked by the reverse shock (corresponding to dou- 
bly shocked external medium originating at r < Rq), and (iv) the 
unperturbed portion of the BM solution which has not yet been 
shocked by the reverse shock (i.e. singly shocked external medium 
originating from r < Ro). These regions are denoted by subscripts 
T through '4', respectively. Regions 2 and 3 are separated by a 
contact discontinuity. 

Immediately after to [or more precisely, at < (f-fo)//o ^ 
a"'^^] the reverse shock reaches only a very small part of the BM 
profile (just behind the contact discontinuity) which corresponds 
to values y - 1 ^ 1 of t he self similar variable, x (defined in 
iBlandford & McKee 1976), so that at this early stage the conditions 
in this region may be approximated as being constant with the val- 
ues just behind the shock for the BM profile with (A,k) = (Ao, Ao) 
at / = /o (i.e. when the shock radius is R = Rq). Since we are in- 
terested in a small range in radius, AR ^ Ro, we can use a pla- 
nar geometry and solve the relevant Riemann problem. Region 
4 is described by the BM solution while in region 1 we have 
pi = wi/c^ = ninip = Air"*', pi = ei = and 71 = 1, where lUp is 
the proton mass. The pressure p, internal energy density e, enthalpy 
density w, rest mass density p, and number density n are measured 
in the proper frame (i.e. the fluid rest frame). We consider a rel- 
ativistic afterglow shock, 74 ^ 1, and a density contrast which is 
not too large such that even after the afterglow shock crosses the 
density bump it will still be relativistic (i.e. 74 ^ t/) ~ a'^*, see 
equation (^J). Under these conditions, in regions 2 and 3 the fluid 
is relativistically hot, p2C^ ^ pi and p^c^ <C P3- Therefore, the adi- 
abatic index in regions 2, 3 and 4 is 4/3, implying p, = e,/3 = Wi/4 
in these regions. This leaves eight unknown quantities: 7, n and 
e in regions 2 and 3, as well as the Lorentz factors of the reverse 
shock, r, , and of the forward shock, F/ . Correspondingly, there are 
eight constraints: three from the shock jump conditions at each of 
the two shocks, and two at the contact discontinuity: 62 = £3 and 
72 = 73. The shock jump conditions simply state the conservation 
of energy, momentum, and particle number across the shock, which 
is equivalent to the continuity of their corresponding fluxes. At the 
rest frame of the shock front they correspond to the continuity of 
wj^ V, W7^ (v /c)^+p, and n'yv, respectively, across the shock (where 
V is the fluid velocity measured in that frame, while p, n, and w are 
measured in the rest frame of the fluid). Unless stated otherwise, 
all velocities and Lorentz factors are measured in the rest frame of 
the unperturbed external medium, which is identified with the lab 
frame where the flow is spherical. 

Under the above assumptions and for 72 = 73 ^ 1 we obtain 



3a-4 



K«-l)-l 



(3) 



In the limit of a relativistic reverse shock (a ^ 1) Eq. j3j reduces 
to V = 74/73 -(30/4)'/". 

The Lorentz factor, 7, of the fluid behind the forward shock 
as a function of ^ R = R/Rq is ^{R < 1) = 74^''"*»>/- before the 



^ For convenience we work throughout the paper in dimensionless vari- 
ables. Unless specified otherwise, x = x/x{Rq). 
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^demonstrates that despite the simphcity of our analytic approxi- 
mation, it provides an excellent description of the accurate solution. 

Once the reverse shock reaches x ^ 2, it samples most of the 
energy in the BM radial profile. At this stage a good fraction of 
the total energy is already in region 2, and a similar energy is in 
region 3. A rough estimate for the radius, /?i = i?o + AA, at which 
this occurs may be obtained by using the conditions in region 2 that 
have been calculated above according to the shock jump conditions 
for a uniform shell, and checking when most of the energy will 
be in region 2. This occurs after a mass ~ iIP'Mq is swept from 
r > Ro, i.e. when the two terms in Eq.|S|become comparable, which 
corresponds to (^'"*' - 1) = a-^ip'Ull -4ki)/(ll -4ko)] ~ a"'^^ or 



^_^tl!' f \1-Aki 



a \ll-Ako 



1/(3-*]) 



(7) 



Figure 1. The Lorentz factor of the forwai'd shock as function of radius 
for three different density profiles that are described by Eq. Q (see leg- 
end for the parameters of each profile). The solid Unes show our analytic 
approximation (Eq. |6l) while the dots are the results of a hydrodynamic 
simulation. 



density jump and 7 = tp^^^^ immediately after the density jump. A 
simple and useful analytic model for 7(.^ > 1) is obtained using the 
energy conservation equation while replacing the mass collected 
up to Ra, Mo = 47rAo/?o ^''/(S-fco), by an effective mass Mcff.o = 
i/j^Mo. The reasoning behind the factor of 1/;^ is to account for the 
fact that just after the density jump the bulk Lorentz factor and the 
average particle random Lorentz factor in region 4 (74) is a factor 
of i/i higher than that in region 2 (72 = 73), so that the energy in 
region 4 is ~ M074 = ip^Moj^. As a result the expression for energy 
conservation ntR>Ro is approximated by 



E = ^Coip^Mo+CiMi (i?)] 7^(i?)c^ = const , 
where Q = 4(3-/t,)/(17-4A:,) (this is valid for fc, < 3), and 



Mi = / 47rr'rfrpex,(/') =Mo f ^— ^ ] a f^"' - 1 



(4) 



3-ki 



(5) 



According to our simple model. 



y(R > 1) = 74 



17-4/to 
17-4/ti 



a[R 



-1/2 



(6) 



Figure Q shows a comparison of our simple analytic expression 
with the results of a hydrodynamic simulation (see appendix lAlfor 
the simulation details). It depicts the forward shock Lorentz factor 
as a function of radius for a spherical ultra-relativistic blast-wave 
that propagates into three different density profiles of the external 
medium that are described by equation Q. Two of the density pro- 
files are uniform both below and above Ro (ko = ki =0) with density 
jumps of a = 10 and a = 100 at /?o. The third density profile presents 
a wind termination shock, for which ko = 2,ki= 0, and a = 4. Figure 



For a > 1 this simplifies to AR/Ro « a"'^^y3(17-4yti)/[2(17- 
4ko)(3-ki)] ~ fl"'^' < 1. Once E/rc' = Cotp^Mo + CMi = M.ff 
becomes comparable to the mass that would have been swept up at 
the same radius if the outer density profile was valid everywhere, 
47rAii?'"*' /(3-fci), the dynamics approach the new BM self similar 
solution for (A.k) = (Ai,ki). By this time Mcff(/?) is dominated by 
the second term in Eq. |4| and therefore the new BM solution is 



approached when R 

R>Rbm = 2''^^-'"\ 



1 becomes comparable to R ' , i.e. when 



2.2 The Resulting Light Curve 

Here the simplified description of the hydrodynamics presented 
above is used in order to obtain a semi-analytic expression for the 
resulting light curve. We obtain explicit expressions for the three 
most relevant power-law segments (PLSs) of the synchrotron spec- 
trum: V < Um < Vc, I'm < v < Vc, and u > max(i',„, i^r), where u,,, is 
the typical synchrotron frequency and is the cooling frequency. 
The first two PLSs appear in the slow cooling regime (v„, < Vc) 
while the last PLS also appears in the fast cooling regime {v,,, > Vc)- 
Two useful time scales for calculating the observed radi- 
ation are the radial time, T, iR) = t — R/c, and the angular time, 
Te{R) = R/2c"y~. The radial time is the arrival time of a photon 
emitted at the shock front at radius R along the line of sight (at 
^ = 0) relative to a photon emitted at / = at = 0. The angular 
time is the arrival time of a photon emitted at the shock front at 
an angle of = 7"' from the line of sight relative to a photon 
emitted at the shock front at the same radius R along the line 
of sight. For convenience we normalize the observed time by 
To = Ti (Ro) = i?o/4(4-fco)c74, f = T /To. In our simple model the 
radial and angular times are given by 



Tr{R > 1) 



Tg(R > 1) = 



1 + 



+ 



^ V 17-4^1 
4-yto\ / 17-4yto 



(4-ko) (R-l) 



, 4-ki J \ \1-Aki 
2(4 -ko)R 



17-4yt, 



(8) 
, (9) 
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while 



Tr{R <l)=R 



Te(R < I) = 2(4 -ko)K 



A-to 



(10) 



Following iNakar & PiranI i2003h . we express the observed 
flux as an integral over the radius R of the forward shock. It is 
convenient to express the integrand, which represents the contribu- 
tion from a given radius R to the observed flux at a given observed 
time T, as the product of two terms^: the total emissivity (per unit 
frequency) of the blast- wave between R and R + dR, A^(R), and 
a weight function, g(T,P), where /3 = dlogF^/dlogv is the spec- 
tral slope, which takes into account the relative contribution from a 
given radius R to the observed flux at a given observed time T: 



F,{T)=Cm 



f 

Jo 



dRA.8(T,f3). 



(11) 



For convenience all variables are normalized by their value at 
To or Ro: K{f) = FAT = fTo)/FATo) and A^R) = A^R = 
RRo)/A^(Ro). The normalization constant C(f3) may be obtained 
by the requirement that Fv{T = 1) = 1, while ^raax(7^) is given by 
Tr[Rmiix(T)] = T and may be obtained by inverting equation |H| for 
R. The weight function g{T,j3) depends on the the dimensionless 
"time" variable 



t{RJ)-. 



T-T,(R) 
fe(R) 



(12) 



and on the PLS where the observed frequency u is in, which is 
specified by the value of the spectral index (3 {F^ oc u^). In princi- 
ple, during the self-similar phase, foi v < Vc, g has a complicated 
form and it also depends on the p ower law index, k, of the exter- 
nal density iNakar & PiranI Eo03h . However, the whole approach 
that is adopted here (of separating the integrand into the product 
of Ai, and a relatively simple g where the time dependence is only 
through r) is anyway not strictly valid during the non-self-similar 
phase. Therefore, we choose to use the simplest expression of g at 
all PLSs, which is the expression obtained in the thin emitting re- 
gion approximation (which is accurate for > Uc), 



AAR<1)=^ p2-3„+*o(3p-5)/4 



^z-.„«op;")/'. <v<Vc, (14) 



For 1 < R < Ri the reverse shock is crossing the hot BM shell 
and is approximated by evaluating the contribution from the 
three emitting regions: A^,2 ~ the shocked external medium origi- 
nating from R> Ro,Av.3 ~ the portion of the BM shell that has been 
shocked by the reverse shock, and - the unperturbed portion 
of the BM solution which has not yet been crossed by the reverse 
shock. 

Regions 2 and 3 can be approximated as being uniform with 
the Lorentz factor given by equation j6). For a relativistic reverse 
shock, which occurs for a large density bump {a ^ 10), we have 



The ratio of the density behind the for- 



,1/2 



tp (3fl/4)'/'' 

ward shock, n2 ~ 473 to that behind the reverse shock, (13 
4(74/273)474^0 = 8(7|/73)no, is n2/'i3 ~ a/ltp^ ^ \/a/3 ~ a' 
Since when both shocks are relativistic the velocity of both shocks 
relative to the contact discontinuity is the same, c/3, the width and 
the volume of regions 2 and 3 is the same (in the contact disconti- 
nuity rest frame), so that their (proper) density ratio is also the ratio 
of the rest mass and the total number of electrons Ne in the two 
shocked regions. Since u,,, oc 767,^, and 7,,, cce/n where €2 = £3 we 



have 7,„2/7m3 =«3/«2 : 



' a and Una/fmj ~ 3/a ~ a 



where B oc e^^e^^^ is the magnetic field measured in the fluid rest 
frame and it is assumed that the fraction of the internal energy in 
the magnetic field, eg = B^/Sire, is the same everywhere. More gen- 
erally, for any value of a > 1, we have 



a(3a + il; ) 



"4 



a 



(15) 



and therefore 



a(a- 



3a + ip^ 



(16) 



g(r,(3) = (l+rf-\ (13) 



An approximation for A„(^) is obtained by considering three 
different phases of emission: ^ < 1, 1 < < J^i, and^ > (similar 
to the approach taken by Pe^er & Wiiers 2006). For ^ < 1 the blast- 
w ave is self-similar and A^iR < I) is similar to the one calculated 
in lNakar & p'irM]<2003h : 



^ This separation is accurate during the self-similar phase (see 
[Nakar & Piran 2003) and serves here as a useful approximation. 



Since_F^,max oc 'yBNe we have F^.max3/-F'„.max2 = 113 /ni = « 

\f3Ja ~ fl"''^^. Therefore, for v < v,„ < Uc or v„, < v < Vc, F^, oc 
Fv.mi.y.Vm and 

d^^(l<^<^0=f-V""=r'/'. (17) 

A^2 V"3/ 

For V > vimx{v,„,Vc) the ratio of the spectral emissivity between 
regions 2 and 3 is equal to the ratio of the energy flux through 
the forward and reverse shocks, respectively, that goes into elec- 
trons with a synchrotron frequency close to the observed frequency, 
i'syn(7e) ~ V- Sincc the magnetic field in the two regions is similar 
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PLS 


M 


f 




m 


7 


P 


V < U,„ < Uc 


1 







-1/3 


2/3 


1/3 


U„, <U <Uc 


1 





(P 


-l)/2 


2p 


(p+l)/4 







2 


(P 


-l)/2 


2p 


(p + 2)/4 



Table 1. The values of the exponents in equations <20l and 1211 for different 
power law segments of the spectrum. 



and so is the total energy flux , the same electron Lorentz factor 7^ 
is required for i^syn(7f) ~ t^, and 



a fluid element within the tail of the BM profile, for which 
'^m{x)/'^m(x = 1) = ;^-<"-5*,)/6(4-t,) jj^^ spcctral emissivity 
per electron scales as P.,.,„^(x)//',...,max(x = 1) = x"<""'*' 
Since during the self-similar evolution x oc i?'^*' while the number 
of emitting electrons in region 3 is constant (Nc, = Mo /nip), we 
obtain 



-7/2 



r'-''-i + (r\-''-i) C-'^Uf- 



(21) 



(1 <R<Ri)-- 



\ P-2 

>..3 \ ^ ^(p-2)/2 

7«..2y 



where 



> max(fm, z^c) . 



(18) 



Region 4 contributes only lo 1/ < Uc- The conditions in this region 
still follow the BM solution, since it does not yet know about the 
density jump, but the fraction of the BM profile that has still not 
passed through the reverse shock decreases with time. This fraction, 
/, is 1 at ^ = 1 and ~ at ^1 . We parameterize f{R) using a linear 
transition with radius. 



(19) 



Summing the contributions from all the different re- 
gions and evaluating the contribution from region 2 as 
Fu oc F^.maxf,'" OC M" 'j'' p^^fR'' we obtain: 



A.(l <R< Ri)g(r,f3) = eiUc-v)f{R)A,,ii^i{R)gi,^Jr,P) 

1 -7/2 

a{R'-'^-\)\ gji^,(T,P), (20) 



17-4*0 
17-4*1 



where 0(x) is the Heaviside step function, and the values of the 
exponents for the relevant PLSs are given in TableQ The subscript 
.R < 1 or > 1 means that the expressions for these R values should 
be used, even if the actual value of R does not fall within this range. 

At the third phase, R > Ri, the reverse shock has finished 
crossing the hot BM shell so that only regions 2 and 3 contribute 
to the emission. Region 2 gradually relaxes into a self-similar 
profile, while region 3 expands and cools at its tail. The con- 
tribution from region 2 remains the same as in the previous 
phase (i.e. at Ro < R < Ri). Region 3 does not contribute at 
V > Vc, while below its contribution can be approximated 
by assuming that its hydrodynamic evolution follows that of 

The energy flux is equal in the limit of a relativistic reverse shock, where 
the velocities of both shocks relative regions 2 and 3 is the same (t/3). 
However, even in the limit of a Newtonian reverse shock (a — I <C 1) the 
velocity of the reverse shock relative to region 3 approaches c/\/3 (the 
sound speed) which is only a factor of \/3 larger than the velocity of the 
forward shock relative to region 2 (c/3). 



A = - 7 [29 - 7yti + ot(37 - 5ki )] 
6 



(22) 



and the power-law indices for the three different PLSs is listed in 
TableQ 



3 CASE STUDIES OF SPHERICALLY SYMMETRIC 
JUMPS IN THE EXTERNAL DENSITY 

In this section we study two spherically symmetric external density 
profiles which are of special interest for GRB afterglows: a wind 
termination shock, and a density jump in a uniform medium. 



3.1 A Wind Termination Shock 

The semi-analytic model for the light curve that has been developed 
in ^2.2l is now applied to a wind termination shock, for which ko = 2, 
ki = 0, and a = 4. We also compare the results of this semi-analytic 
model to the numerical model that is described in Appendix IaT. 
The resulting light curves are displayed in Figures l2m for the three 
most relevant power law segments (PLSs) of the spectrum. The re- 
sults of the semi-analytic model nicely agree with the numerical 
results. Some differences do exist but the qualitative behavior (i.e. 
variation time scales and amplitudes) is similar, and even the quan- 
titative differences are not very large. The main differences between 
the semi-analytic model and the numerical simulations are a small 
initial dip before the rise in the flux for f < f^, a difference in 
the exact starting time for the change in the temporal decay index 
for Urn <u <Vc and v > m&TL{Um,Vc), and a slightly different nor- 
malization of the asymptotic flux at T ^ To for u <i/c- The latter 
arises since we neglect the dependence of the function g on A; in 
these PLSs (g does not depend on kfoTv> Vc) in the semi-analytic 
model. This causes a deviation (by a factor of the order of unity) 
in the normalization of the asymptotic flux calculated by the semi- 
analytic model, compared to its true value, in cases where ko=^k\. 

The light-curves show a smooth transition between the asymp- 
totic power law behavior at T < To and at T ^ To. There is no 



This code assumes optically thin synchrotron emission and does not take 
into account opacity related effects such as synchrotron self-absorption or 
synchrotron self-Compton. 
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Figure 2. The light curve {upper panel) and the evolution of the temporal 
decay index a (lower panel) for the synchrotron emission from a spherical 
relativistic blast wave running into a wind termination shock (Uq = 2,k\= 0, 
fl = 4 in Eq. (3), in the frequency range v < v,n < Vc- Shown are the results 
of the semi-analytic model described is ^2.2V solid thick blue line) and of the 
numerical code described in AppendixlAlffc/ac^ dots) for a wind termination 
shock. The semi-analytic result for the same wind without a termination 
shock (where pcxt = Aof"*" at all radii; red thin line) is added for reference. 



rebrightening in PLSs where the flux decays at 7 < To (i^ > Vm), 
and no sharp feature in the light curve which might serve as a clear 
observational signature. The transition between the two asymptotic 
curves {T < To & T ^ To) is continuous with the temporal decay 
index rising slowly for u <Vc and mildly fluctuating for > Vc. The 
values of the temporal index a = d\o^F,j / dlogv during its rising 
or fluctuating phase do not exceed its asymptotic value at T 3> To by 
more than 0.1, at any time. Therefore, the only observable signature 
of a wind termination shock is a continuous break (with Aa = 0.5) 
below i^c- Above Vc there is no change in the asymptotic value of 
the temporal index a, and it only slowly fluctuates with a very 
small amplitude (~ 0.1), which is extremely hard to detect. Note 
that these results are applicable for a case where the blast-wave re- 
mains relativistic also after it encounters the termination shock (i.e. 
73 ^ 0.7274 > 3). 

The break in the light curve (the shallowing of the flux decay 
for v,„ < V < Vc OT the transition from constant to rising flux for 
< < Vc) occurs over about one decade in time. Initially there 
is very little difference relative to the case where there is no wind 
termination shock (and p^xi = Aor"*" at all radii), or even a small 
dip ally > max(iy„,iyc), while a rise in the relative flux starts at T = 
T /To ~ 2-4. The light curve approaches its asymptotic late time 
power law behavior at about f ~ 10- 10^. This can be understood 
as follows. 

The contribution to the observed flux from within a given an- 
gle 9 around the line of sight does not change drastically across the 
density jump. However, there is a sudden decrease in the Lorentz 
factor of the shocked material behind the forward shock, so that the 



Wind termination sliock 
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T/To 

Figure 3. Same as Figure l2lbut for the spectral range <v < Vc- We use 
p = 2.5. 
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Figure 4. Same as Figurel3lbut for the spectral range u > max(Ui„,Uc), and 
with an additional panel for the ratio of the flux with and without a wind 
termination shock. 
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effective visible region increases from 6 < 74' io 9 < 73' = ip/'yA- 
This is responsible for most of the observable signature, and it starts 
affecting the light curve noticeably only when photons emitted just 
after Ra from an angle 9 > 74' arrive to the observer, namely at 
~ TeiRo.-) = 4To where 



To(Ro.±) ■: 



\imTel{l±e)Ro] 



(23) 



This full angular effect becomes apparent when photons from 
the same radius and an angle of ~ 73' reach the observer, 
at ~ Tg(Ro.+) = 2(4 - ko)il>^To ~ STo. The radial time is smaller 
than the angular time, and therefore the radial effect would only 
slightly increase the time when the total effect becomes promi- 
nent. The light curve approaches its asymptotic power law be- 
havior when the dynamics approach the new self similar evolu- 
tion, at ~ ^BM = 2'^*^"'''. Since [74/7(^81^1)]^ ~ 4 and therefore 
7e(^BM) ~ 16, while TriRau) ~ 2.6, this corresponds to T ~ 20. 
Obviously, this is a rough estimate, but it agrees reasonably well 
with the numerical results. 

Our main result for the light curves from a wind termination 
shock is that there is no prominent readily detectable signature in 
the light curve. This is very different from th e results of previous 
papers that explored a wind termination shock ( Ramirez-Ruiz et^lj 
2001, 2005; Dai & Lu 2002; Eldridge et al. 2006; Pe'er & Wij^ 
2006) which predicted a clear observational signature (including 
optical rebrightening also when the termination shock is at a suffi- 
ciently small radius so that the blast-wa ve is still relativistic when 
it runs into it). In some of these works ( Ramirez-R uiz et alj200lL 
[2005) the forward shock becomes non-relativistic after hitting the 
density jump, which might account for some differences compared 
to our results which are valid for the case when the forward shock 
remain s relativistic after running into the density jump. In othe r 
works jPai & L u 2002 ; Eldridge et al. 2006; Pe'er & Wiiers 200^, 
however, the forward shock is assumed to remain relativistic after 
encountering the density jump, similar to our assumption. The main 
reason for the discrepancy relative to the latter works is that they did 
not consider either t he effect of the reverse shock on the dynamics 
jEldridge etalT iooe) or did not properly account for the effect of 
different photon arrival times fro m different angles relative to the 
line o f sight from any given radius jPai & Lj200ilPe'er & wiier^ 
l2006h . Both effects tend to smooth out the resulting variability in 
the light curve. 



3.2 Spherical Jump in a Uniform External Medium 

Next we explore the light curve that results from a spherical rela- 
tivistic blast-wave running into a uniform external density with a 
jump at some radius Ro (i.e., ko =ki =0 and a > 1). Such a density 
profile can be generated, for example, by the contact discontinuity 
between shocked winds from two evolutionary phases of the mas- 
sive star progenitor. This configuration also serves as an approxi- 
mation for a large density clump, and constrains the observational 
signature from a small density clump (see 

Figures 15 17 1 depict the light curves from our semi-analytic 
model (described in ^2.2> for four different density contrasts (a = 
2, 10, 100, 1000), as well as the results of the numerical simulation 
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1^ 




~^^T^^^ ' ' ' ' 
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T/T, 



10 



10 



Figure 5. The light curve (upper panel), the ratio of the flux with and with- 
out a density jump (middle panel), and the evolution of the temporal index 
a (lower panel) for the synchrotron emission from a spherical relativis- 
tic blast wave propagating into a medium with a step like density profile 
(Icq = k\ = in Eq. Q), in the frequency range u < u„, < v^. Shown are the 
results of the semi-analytic model described is i|2.2l for four different den- 
sity contrasts (fx = 2, 10, 100 & 1000) and of the numerical code described 
in AppendixlAlfor two cases (a = 10 & 100). We use here p = 2.5. 



(described in Appendix]^ for two of these cases (a = 10, 100). 
The agreement between the semi-analytic model and the results of 
the simulation is satisfactory. In all cases (all the different PLSs 
and a values) the semi-analytic model qualitatively follows the nu- 
merical results and recovers the main features (i.e. the correct time 
scales and amplitudes of the variations and their derivatives). In 
most cases the quantitative comparison is also impressive (better 
than 10%). The main differences between the semi-analytic model 
and the simulation results are the small initial dip before the rise 
in the flux for u < v,„ < i^c (observed also in the wind-termination 
sock) and an over-shoot for a = 100 in this PLS. The semi-analytic 
model predicts an initial dip for u,„ <v<h'c, which also appears, 
although less prominently, in the results of the numerical simula- 
tions. 

The main result that emerges from Figures I5I/I is that no 
sharp features appear in any of the light curves, no matter how 
high the density contrast (at least as long as 73 = 74/1/) 3> 1 where 
-i/j ~ (3a/4)'^'' for a ^ I). Moreover, the maximal deviation of the 
temporal decay index, q, from its asymptotic value (which is the 
same for T < To and T 2> To, since = ki) is not large (< I in all 
PLSs at all times), and as we show for iy,„ < v < Pc it approaches 
an asymptotic value at large a. 

Observationally, the most interesting PLS is v,„ < v < Vc, 
since it typically includes the optical. In this PLS the flux enhance- 
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Figure 6. Same as Figure|5]but for the spectral range v,n < u < Uc. 



ment, f(T) = FAT)/FATm = 1), is asymptotically f(T > To) = 

1 /2 

a ' , but the transition to this asymptotic value is very gradual. 
Figure |8| depicts the value of the maximal deviation of the tem- 
poral index a from its asymptotic value, Aomax, as a function of 
a for two values of p. We emphasize the behavior of Aa since it 
is perhaps the easiest quantity to observe. The inability of density 
fluctuations to produce sharp features in the light curve is demon- 
strated by the low values of Aoniax that we find. Some examples are 
0.1, 0.4 (0.35), 0.75 (0.65) & 0.95 for a = 2, 10, 100 & 1000, re- 
spectively, where the values in brackets correspond to the results of 
the numerical simulation. Furthermore, for very large values of a, 
Aomax saturates at a value of ~ 1. Since a{T < To) = -3(p- 1)/4 ~ 
-1, no rebrightening (i.e., a > 0) is observed. Moreover, at first 
(just after To) a mild dip, is apparent in the light-curve. The depth 
of this dip increases with a. Another constraining observable is the 
time over which Aa„,ax is obtained. Thus, we consider the ratio of 
the time when Aomax is obtained and the time when Aa > once 
it recovers from its initial dip. We find this time ratio to be « 5 for 
any reasonable values of a and p. 

Our results can be understood as follows. The contribution 
to the emission per unit area of the shock front along the line of 
sight, from a radius R to an observer time T, increases with the 
energy density of the shocked fluid and with its bulk Lorentz fac- 
tor. A density jump immediately increases the energy density of 
the freshly shocked fluid (by a factor atp~^ > 1) while reducing its 
Lorentz factor (by a factor t/j^^ < 1). The net effect is such that 




Figure 7. Same as Figure l^ut for the spectral range Vm^Vc <v. For clarity 
we omitted the a = 2 curve. 

the decrement in the Lorentz factor dominates by a small margin, 
and the line-of-sight emission actually drops at Rq. This drop in- 
creases with a and is the source of the observed dip right after To. 
The same drop in 7 is also the origin of the flux increase that fol- 
lows. With a lower 7 the largest angle 6 (from the line of sight) 
that contributes to the observed emission increases. This contri- 
bution becomes apparent only when emission from 6 > 1 /74 ar- 
rives to the observer, at ~ Te(Ro.-) = STo, and is completed when 
emission from 6 > I /^j, is observed, at ~ Tg{Ro.+) = ip'TeiRo-) ~ 
o^^'TgiRo-). The transition to the asymptotic value continues up 
to Tbm ^ a''^^Tg(Ro.+) ^ lOaTo. These time scales explain why 
the transition is so gradual, and so is the convergence of a to its 
asymptotic value at T 3> To, for large a. Our results show that 
the maximal value of a is observed around Tg(Ro.+), and that for 
fl > 1, flTg(Ro.+)] ~ (0.1-0.4)fl'/l Now Aa can be approximated 
by log{f[Tg{Ro.+)]}/log[Tg(Ro,+)/Tg{Ro.-)] which approaches 1 at 
large a. 

It is generally accepted that the main signature of density fluc- 
tuations in the external medium are chromatic fluctuations in the 
afterglow light curve, where sharp features are expected below i^c 
(which typically includes the optical bands) and no (or very weak) 
variability is expected above (which typically includes the X- 
rays). This conclusion relies on the fact that a change in the exter- 
nal density effects the asymptotic light-curve (at T ^To compared 
to r < To) only below Vc, but not above Uc- A comparison between 
Figures|6|&Qshows that while the behavior in these two PLSs is 
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Figure 8. The maximal deviation of tlie temporal index (a = 
d\ogFv / dlogT) from its asymptotic value, Aomax, as a function of the 
density contrast a, for electron power-law indexes p = 2.5 (solid line) 
and p = 2.\ (dashed line). The results of the numerical simulations for 
a = 10 & 100 (in which p = 2.5) are marked as dots. 



indeed different, the general concept described above is inaccurate 
and the differences are more subtle. Both PLSs show smooth fluc- 
tuations in the temporal index a with a comparable amplitude. The 
main difference is in the flux normalization of the asymptotic light 
curve at r 3> 7b compared to T < To- Above Vc the asymptotic light 
cvu^e does not change and the observed flux fluctuates around this 
asymptotic power law decay, while below Vc the normalization for 
the asymptotic light curve at T 3> To is larger than that at T < To by 
a factor of a'''", and therefore the flux continuously increase rela- 
tive to the case where there is no density jump (a= 1). This type of 
difference in the behavior is, however, much harder to detect (com- 
pared to variations in a) since, obviously, the reference light-curve 
(for a=i) cannot be observed. 



ing from such a density profile I 


Lazzati et alJ l200l iDai & Ld 




iNakar. Piran & GranotlflOOS 


Nakar & Pirarl 


l2003h and all 



of them predicted much sharper features with an observable re- 
brightening at I'm < V < Vc (i.e. a > 0) which doe s not exist in 
the results presented here. In several wo rks I Lazz atTet alJl200^ 
iNakar. Piran & Granotll2003l:lNakar & Piraa.20 03) the main cause 
for the discrepancy is that the effect of the reverse shock on the dy- 
namics was neglected. As we show here, even if the reverse shock 
in only mildly relativistic (a ~ 2), its effect on the dynamics cannot 
be neglected. The reason for this is that even if the emission from 
the reverse shock itself is negligible, the abrupt drop in 7, which is 
the Lorentz factor of the shocked material behind both the reverse 
and the forward shocks, prevents ver y signific ant variations in a 
also from the forward shock emission. lPai & Lu l 2QQ2) included a 
partial consideration of the reverse shock, however they neglected 
the strong effect that angular smoothing has on the light curve. 



3.3 The Effects of Proximity to a Break Frequency 

So far we have assumed that the observed frequency 1/ is very far 
from the break frequencies (v,,, and Vc) and therefore remains in 
the same power law segment (PLS) of the synchrotron spectrum 
throughout the hydrodynamic transitions that we have investigated. 
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Figure 9. The temporal evolution of the typical synchrotron frequency Vm 
(dotted line), which is defined as the frequency at which the spectral in- 
dex 13 is midway between its value at <C v„, (I3[ = 1/3) and at f 2> i^m 
[P2 = (1 -p)/2], /3 = (/3i +/32)/2. We use p = 2.5. The shaded region shows 
the frequency range where 80% of the change in (3 occurs (between 10% 
and 90% of A/3 = /3i —fii). The frequency where the asymptotic power laws 
well above and below v,,, meet is marked with filled squares. For clarity, all 
frequencies are normalized by v„, q = u„,(Tq) and multiplied by the appro- 
priate power of T = T /To which takes out the asymptotic time dependence 
of Um at early and late times. The different panels are for different spherical 
density profiles with a sharp density jump that have been studied in il3.1l - 
a wind termination shock (upper panel), and in il3.2l - a spherical density 
jump by a factor of a = 10 (middle panel) and a = 100 (lower panel) in a 
uniform medium. 



Under that assumption the flux density normalized by its value at 
To is independent of frequency within each PLS. In this subsection 
we examine the effect on the light curve if the observed frequency 
is in the vicinity of a break frequency around the time of the hydro- 
dynamic transition. For this purpose we use our numerical code* 
and consider the cases that have been studied numerically in ^3.11 
(a wind termination shock) and in ^3.2l (a spherical density jump 
by a factor of a = 10 or a = 100 in a uniform medium). 

Figures |9| and 1101 show the temporal evolution of the spec- 
tral break frequencies v,,, and Uc, respectively, around the time 
of the density jump. The break frequencies are defined as^ i/j = 



* Since our semi-analytic model was designed only for the case where the 
observed frequency remains in the same PLS, it is not appropriate for this 
purpose. 

^ Defining instead Vj, = max{j^|/3(j^) > 0.5/3i +0.5/32} makes no noticeable 
difference. 
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Figure 10. The temporal evolution of the cooling break frequency, Uc, in the 
same format as Figure l9l here (3\ = (i-p)/2, Pi = -p/2, and again p = 2.5. 
For a wind termination shock {upper panel) the asymptotic temporal index 
at T < To (qi = 1/2; left dashed line) is different from that at T 2> To 
(a2 = right dashed line), and therefore the normalized frequencies 

are not multiplied by (T /Tq)^!^ as in the other two panels (that are for a 
density jump in a uniform medium for which cji =012 = -1/2). 
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Figure 11. The light curves for frequencies that are in the vicinity of Um at 
the time of the density jump, for the same hydrodynamic transitions that are 
shown in Figure|9](with p = 2.5). The solid thick line is for u i/,„(7o), the 
dashed thick line is for v ^m.lo%{To)' and the doited thick line is for v 
^m.90%(To)- The remaining curves are spaced by a factor of 2.5 in frequency 
in the middle panel and 3 in the top and bottom panels. 




min{z/|/3(z/) < 0.5/3i + 0.5/32}, where /3i and (32 are the asymp- 
totic values of the spectral index Patv<^ vi, and vt, respec- 
tively, and b = m,c. The shaded region shows the frequency range 
1^10% < y < 1^90% where i^io% = min{i^|/3(i^) < 0.9/3i +0.1/32} and 
v<xi% = m?ix{v\l3(v) > 0.1/3i +0.9/32}, i.e. where 80% of the change 
in 13 across the spectral break occurs. 

The typical synchrotron frequency fluctuates around its 
asymptotic T^^^^ power law decay, which does not depend on the 
power law index k of the external density. Furthermore, even for a 
wind termination shock there is no observable change in the asymp- 
totic normalization (i.e. the asymptotic value of T^^'Vm) and only 
a very small change in the width of the spectral break which is de- 
picted by the shaded region (in agreement with the semi-analytic 
results of iGranot & Saril l2002h . For a density jump in a uniform 
medium there is no change in the asymptotic normalization or in 
the asymptotic shape of the spectral break, again as expected from 
analytic calculations. Both the amplitude and the typical time scale 
of the fluctuations in T^^^Vm increase with the density contrast a, 
where the amplitude scales roughly as a'^^ and the time scale is 
roughly linear in a. There is a sharp feature in which occurs first 
for Vm.wk, then for Vm, and finally for f,„.9o%- This can be under- 
stood as follows. Immediately after the forward shock encounters 
the density jump v,„ decreases in region 2 and increases in region 3, 



by a factor of ~ a'l^ in both cases. Therefore, as long as the region 
3 contributes significantly to the observed spectrum, the break is a 
superposition of two peaks (corresponding to t',,,.2 and i/m.s), sepa- 
rated by a factor of ~ a in frequency, and thus its width increases 
significantly. Moreover, it causes ii(v) to be non-monotonic, and 
the definitions of f,„.io%, v,,,, and Vm.wh cause these frequencies to 
have a finite jump when the extremum in /3(i') crosses the appropri- 
ate value of /3, which occurs later for higher frequencies. Given the 
complex structure of the spectral break around f,,, we also present 
in figure |9| the evolution of the frequency where the asymptotic 
power laws well above and below the break meet (which can serve 
as an alternative defi nition for the lo cation of the break frequency, 
as was done in Gran ot & Sarl2002h . This frequency evolves very 
smoothly and shows very mild fluctuations in all cases, since it is 
less effected by the transient broadening of the spectral break dur- 
ing the hydrodynamic transition. 

The evolution of the cooling break frequency, v^, is shown 
in Figure 1101 For a wind termination shock the temporal index 
a. has different asymptotic values at T < To (ai = 1/2) and at 
T ^ To (a2 = -1/2); i/c transitions rather smoothly between these 
two asymptotic limits, with a slight overshoot (i.e. dlogUc/dlogT 
dips below -1/2) due to the increase in density across the jump 
(the asymptotic value of i/c decreases with increasing density). 
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The asymptotic behavior of Vc is marked with dashed lines show- 
ing that for practical purposes (e.g., analytic calculation) it can be 
well approximated as a sharp temporal transition between al and 
02 at r ~ 1.5 To, while keeping in mind that the spectral break 
itself is very smooth at any time. For a density jump in a uni- 
form medium the asymptotic value of the temporal index does 
not change (ai = 02 = -1/2), but the asymptotic normalization of 
r'^^j^r at r ^ To is a factor of a lower than at T < Tq. The transi- 
tion between the two asymptotic limits is fairly smooth. The shaded 
region which corresponds to 80% of the change in [3 across the 
break is significantly larger for i^c compared to v,,,, corresponding 
to a smoother break , in agreement with semi-analytic calculations 
iOranot & Sari 2002). The transition to the late time asymptotic be- 
havior stretches over a larger factor in time, which increases with 
a. 

Figures fTTI and fT2l show light curves for frequencies that are 
in the proximity of v,,, and Vc, respectively, around the time of the 
density jump. There is a smooth transition between the light curves 
for frequencies that are well below the break frequency and those 
for frequencies well above the break frequency around the time of 
the density jump. Our main conclusions from sections § 93. H and 
13. 21 remain valid also when the observed frequency is near a break 
frequency around the time of the density jump. In particular, there 
is no rebrightening ai v ^ v,„, and the observed features in the 
light curve are very smooth. Therefore, our main results are rather 
robust. 



4 A CLUMP IN THE EXTERNAL DENSITY 

In this section we estimate the effect of a clump in the external 
density on the light curve. By a clump we refer to a well localized 
region of typical size /d which is overdense by a factor of a > 1 
relative to the uniform background external density. For a given 
clump size la and overdensity a (well within the clump, near its 
center) the effect on the light curve is expected to be larger if the 
clump has sharper edges, i.e. the smaller the length scale A/ over 
which the density rises by a factor of ~ a relative to the background 
density. While in practice one might, in many cases, expect A/ ~ 
Zci, we consider the limit of a sharp edged clump with A/ ^ la, in 
order to maximize the effect on the light curve. 

Because of relativistic beaming, most of the contribution to 
the observed light curve from a given radius R is from within an 
angle of 6 < 7"' around the line of sight, which corresponds to a 
lateral size of ~ R/j. Therefore, a clump in the external density of 
size^ /ci ^ R/j would not differ considerably from a spherically 
symmetric density jump that was considered in ^ If the surface 
of the clump is not normal to the line of sight (e.g. if the line of 
sight to the central source does not pass through the center of a 
spherical clump), this is expected to reduce the effect of the clump 
on the light curve (similar to what is expected if the clump does not 
have very sharp edges, Al ~ la rather than <C la)- Therefore the 
results of ^can be viewed as a rough upper limit on the effect of 

^ Here 7 is the Lorentz factor inside the clump, which is smaller than that 
before the afterglow shock hits the clump, by a factor ip that is given in 
equation j3). 




Figure 12. The light curves for frequencies near Uc at the time of the density 
jump, for the same hydrodynamic transitions that are shown in Figure [Tol 
(with p = 2.5). The noi-malized flux is multiplied by (T /To)'-^''-^^/* f "5 
in order to eliminate asymptotic time dependence at 7 S> To for a uniform 
medium. The solid thick line is for v i'c(7b), the dashed thick line is for 
u lO%(^o)> th^ dotted thick line is for v ^ t^cmATo)- The remain- 
ing curves ai'e spaced by a factor of 15 in frequency in the top panel and 
8 in the middle and bottom panels. 



"big" clumps {la > Ro/"/3)- Small to intermediate clumps, of size 
Id ^ ^0 /73 are expected to have a smaller effect on the light curve 
and are investigated below. 

The results of the previous section, and in particular the semi- 
analytic model for the light curve that was derived in ^2.21 can be 
used to put an approximate upper limit on the effect that a den- 
sity clump could have on the observed light curve. Such a limit 
is achieved by using the spherical model from 92.21 within a fi- 
nite solid angle: Smin < < O^ax and 0min < 4> < 4>max in spheri- 
cal coordinates, while in the radial direction the clump extends out 
to R ^ Rq. In practice we expect the radial extent of the clump 
to be AR ~ la, similar to its extent in the 9 direction (~ R0A6 
where A6 = Smax - dmin) and in the <j) direction (~ RoOAcp where 
max 0min and 9 = [0min + ^max]/2). In our coarse approxi- 
mation the clump has no upper bound in the radial direction, and 
its lateral size scales linearly with radius. Obviously, this sets an 
(approximate) upper limit for the effect on the light curve of a 
clump with a size la in all directions (which is roughly given by 
la ~ RoA9 ~ Ro9A(j) ~ AR in spherical coordinates). 

Figure fTsl shows the results of this model for a density clump 
that lies along the line of sight, 6'niin = and 74&max = 1/3, 1, 00 
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Figure 13. The ratio of the flux with and without a density jump at S > Rq, 
within a finite angle 9 < Smax around the line of sight, for three differ- 
ent density contrasts (a = 10, 100, 1000) and three different angular sizes 
{'y^dmax = li '^)- This serves as an approximate upper Hmit for the 
effect of a clump in the external medium on the Ught curve. The differ- 
ent panels are for the most relevant power law segments of the synchrotron 
spectrum. 



while A(j> = 2tt, for three different values of the density contrast 
(a = 10, 100, 1000) and for the three most relevant power law seg- 
ments of the spectrum (that were modeled in 82. 2> . The smaller 
the angular extent of the clump, the smaller the amplitude of the 
change in the flux relative to its value for the smooth underlying 
density distribution without the clump (a= 1), and the smaller the 
factor in time over which it effects the flux significantly (e.g., the 
full width at half maximum of the relative flux). This behavior is 
expected since both the amplitude and the duration of the fluctu- 
ation depends on the size of the clump ( loka, Koba vashi & Zhan3 
llOOS *). The amplitude depend on the ratio in size between the per- 
turbed region (of length scale 1^) and the unperturbed region (of 
length scale R/y) of the blast wave, and thus increase with /(■;. The 
duration depends on the delay in the arrival time of photons emit- 
ted from the perturbed region, which again increase with the size 
of this region. 

As can be seen in Figure fTsl the amplitude of the fluctuation 
in the relative flux increases with the density contrast a. There is a 
sharp transition in the light curve at J" = 1 -H2(4-A'o)(746max)^, when 
the radiation from the outer edge of the clump at R = Rq reaches 
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Figure 14. Similar to Figure lTsl but for a fixed density contrast, a = 100, and 
for three clumps along the line of sight together with a clump that is at the 
side of the visible region at the time of collision, centered around 746 Ri 1. 



the observer, which corresponds to the peak in the relative flux for 

= 746'maxi/' ' ~ 1- Such sharp features (see also Figure 
are caused by the over-simplified clump model that we use here, 
and are expected to be smoothed out in more realistic models of 
clumps. Above i/,,, a clump can produce a dip or a bump in the 
relative flux with a relatively small amplitude (depending on its 
size and density contrast), while below 1/,,, it produces a bump in 
the relative flux with a larger amplitude. 

Figure fT4l shows the relative flux for a fixed density contrast, 
a = 100, for three clumps along the line of sight with different an- 
gular sizes (A</i = 27r, Smin = and 746inax =0.1, 1/3, 1), as well 
as for a clump close to the edge of the visible region at the time 
of the collision, 74(0min, ^max) = (2/3,4/3) and Acp = n/6, which 
occupies the same solid angle as the clump along the line of sight 
with 74Sraax = 1 /3. The effect of a given clump on the light curve is 
maximal when it is along the line of sight, and it becomes smaller 
the further it is from the line of sight (i.e., its effect is significant 
over a shorter time scale and the amplitude of the change in the 
relative flux is somewhat reduced). 

Overall, a fairly large clump of size la ^ R/'y with a suffi- 
ciently large density contrast (a > 10- 10^) is required in order to 
produce an observable signature in the light curve, and even then 
the most prominent signal would be below i/,,,, which is typically 
relevant for the radio. In the optical, which is typically above Um, 
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there would only be a small dip or bump in the relative flux, which 
would be hard to detect. 



5 IMPLICATIONS TO GRBS 021004, 000301C, 030329 
AND SHB 060313 

Next we reconsider the cause for the fluctuations in the optical af- 
terglow light curves of the long-soft GRBs 021004, 000301C, and 
030329, as well as the very recent short-hard GRB 060313, in light 
of our results. GRB 021004 showed significant variability in its op- 
tical afterglow (Pandevetal. 2002; Fox et al. 2003; Bersieretal. 
I2OO3 ': Uemura et al. 2003) which included three distinct episodes 
of mild rebrightening (a > 0) at T ~ 0.05 days, ~ 0.8 days, and 
~ 2.6 days. Just before the first of these epochs, a ~ -0.7, and it 
then became positive over a factor of ~ 2-3 in time. Such a large 
increase in a over a relatively small factor in time is hard to ac- 
commodate by a jump in the external density: Aa ~ 1 requires 
a > 10^ and even then it is hard to achieve such an increase in a 
over a factor of ~ 2 - 3 in time. The second rebrightening episode 
lies within the tail of the first, and is therefore not a very "clean" 
case to study. The third rebrightening epoch at ~ 2.6 days is some- 
what more isolated, and has Aa > 1 where most of the increase in 
Q is within about a factor of ~ 1.5 in time. This is extremely hard 
to achieve by variations in the external density. On the other hand, 
angular fluctuations in the energy per solid angle within the jet (a 
"patchy shell") nicely account for both the fluctuations in the light 
curve and the variability in the linear polariza tion of this afterglow 
jGranot & K6nigl2003l:lNakar & Orenl2004h . 

GRB 00030 IC display ed a largely achromatic bump in its 
optical to NIR light curves ('Sagar et al.'' 2000l : iBerger et alJ l2000') 
peaking at T ~ 3.8 days. Before the bump a « -(1.2- 1.3) and 
during the rise to the bump a became slightly positive, so that 
Aa ~ 1.5, where most of the increase in a occurred over a fac- 
tor of 1.5 in time. We find that this cannot be produced by 
a sudden change in t he external density (as has been suggested 
bv lBergeretaljboOd) . The decay just after the peak of the bump 
is ver y sharp, a « -(3.5-4), and since /3-q > 2 iSagar et'all 
I2OOOI find P = -0.96 ± 0.08 during the decay just after the peak, 
at T = 4.8 days) this r equires a deviation from spherical symme- 
try teumar & Panait escu 2000), and might not be easy to achieve 
with angular fluctuations in the energy per solid angle within the 
jet (a "patchy shell") or aspherical refre shed shocks. In this case an 
altern ative microlensing interpretation taamavich. Loeb & Stanekl 
I2OOOI) was shown to be a ble to nicely reproduce the shape of the 
bump jGaudi. Granot & L oeb 2001). 

GRB 030329 has one of the best monitored and densely sam- 
pled opti cal afterglow light curves to date. It is presented in great 
detail bv iLipkin et alj i2004h who clearly show that the light curve 
contains several rebrightening (a > 0) episodes in which Aa > 2. 
All of these episodes have a similar structure and occur over a small 
factor in time (AT < T). In previous works it has been suggested 
that some o f these rebrightening ep isodes are a result of density 
fluctuations. IPe'er & WiiersI yP06l have suggested that the first 
(and largest) rebrightening episode is the signa ture of a wind termi- 
nation shock. Sheth et al. 1 2003|) have foUowed lierger et aljiioosi) 
in attributing the first rebrightening episode to a two-component 
jet, but argued that the subsequent rebrightening episodes could 



be explained by variations in the external density. Our results 
clearly show that none of the bumps in the optical light curve of 
GRB 030329 can be a result of density bumps. This implies that 
the two-jet model of Berger et al. (2003) and Sheth etal. l20o3) 
fails to account for most of the observed light-curve fluctuations. 
Excluding density fluctuations as the possible source of any of 
the maj or bumps, together with the similarity in the shapes of the 
bumps ( Liokin et al. 2004), support a sequence of similar episodes 
in which the energy of the jet fl uctuates, such as a late time en ergy 
injection by "refreshed shocks' ' jGranot. Nakar & Piraj2003l) . 

The afterglow of the very recent short-hard GRB 060313 
has been monitored both by the X-ray telescope (XRT) and 
by the opt ical/ultra-violate telescope (UVOT) on board Swift 
teoming et al...2006.) . The optical/UV light curve showed three 
sharp bumps/flares at T ^ 1.7 hr, ^ 3.2 hr, and ~ 6.6 hr, with an 
amplitude of more than a factor of 2 in flux and a very short rise 
time of Ar < O.IT. During the same time the X-ray light curve 
showed a smooth (and steeper ) power law decay. This was inter- 
preted bv lRoming et alji2006h as the result of variations in the ex- 
ternal density by a factor of ~ 2, where the lack of variability in 
the X-rays was attributed to Vc being between the optical/UV and 
the X-rays. Our results clearly show that such sharp rebrightening 
episodes as were seen in the optical/UV afterglow light curve of 
GRB 060313 cannot be the result of variations in the external den- 
sity. Therefore, there must be some other cause for this variability, 
such as late time internal shocks with a very soft spectrum, as was 
suggested bv Roming et al. (.20061) as an alternative mechanism. 



6 CONCLUSIONS 

We have presented a semi-analytic model for the light-curve result- 
ing from synchrotron emission by a spherical relativistic blast-wave 
that propagates into a power law external density profile with a sin- 
gle sharp density jump at some radius Rq. Our solution is general 
enough to include a transition in the density power-law index {k) 
at ^0, but is limited to cases in which the blast-wave remains rel- 
ativistic after it encounters the density jump. This model has been 
used to explore in detail two density profiles that are most rele- 
vant to GRB afterglows: a wind termination shock, and a sharp 
density jump between two regions of uniform density. The latter 
results are also used to constrain the signature of density clumps 
in a uniform medium. We have also carried out detailed numerical 
simulations for several of the cases which we have studied in detail. 
These numerical results serve three purposes. First, they are used 
in order to obtain more accurate light curves for the cases which 
are of special interest. Second, they are appropriate for calculating 
the light curves in the vicinity of a spectral break frequency ( 33.31 . 
Third, they serve as a test for the quality of our simple semi-analytic 
model, which is found to give a very good qualitative description 
and reasonable quantitative description of the light curve. 

Our main result is that density jumps do not produce sharp 
features in the light curve, regardless of their density con- 
trast! The results of our specific case studies are as follows. 

A wind termination shock: 

• The light curve shows a smooth transition, which lasts for 
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about one decade in time, between the asymptotic power-law be- 
havior at r < To and at T ^ To. 

• There is no rebrightening or any other sharp feature that can 
be used as a clear observational signature. 

• Above the cooling frequency, Uc, there is no change in the 
asymptotic value of the temporal index, a, and it only fluctuates 
with a small amplitude (Aa « 0.1). 

• The only observable signatures of a wind termination shock 
are a smooth break, with an increase of Aa = 0.5 in a, below Uc 
and a transition in the temporal evolution of Vc- 

A density jump between two uniform density regions: 

• The light curve shows a smooth transition between the two 
asymptotic power laws (at T <Tq and T 2> To). 

• The transition time increases with the density contrast, a, and 
is about ~ 10a To. 

• The maximal deviation, Aofmax, of the temporal index a from 
its asymptotic value (at T <To and T 3> 7b) is small. For example, 
AQmax(a = 10) < 0.4; Aamax depends weakly on a and approaches 
~ 1 at very large a values. Therefore, a density jump cannot pro- 
duce an optical rebrightening when /^optical > i-'m- 

• The light curve fluctuates also above (typically including 
the X-ray band). While the asymptotic flux (at T 2> To) above Uc is 
unaffected by the density jump, the fluctuations in a are compara- 
ble to those below Vc- 

An overdense clump on top of a uniform density background: 

• Only a fairly large clump (la > R/'y) with a sufficiently large 
density contrast {a> 10-10^) produces a significant fluctuation in 
the light curve. 

• The effect of a clump on the light curve is significantly larger 
when it is located along the line of sight, than at an angle of ~ I/7 
from the line of sight. 

• The signature of a clump is most apparent al u < < Vc- 

• Above I'll, a clump can actually cause a small dip in the light 
curve, while below i^,„ it causes a (larger) bump. 

For a spherical density jump our conclusions are based on accu- 
rate results, while in the case of a density clump we obtain only 
an approximate upper limit for its effect on the light-curve. There- 
fore, our results for density clumps should be taken only as rough 
guidelines. Our main results remain valid also when the observed 
frequency is close to a spectral break frequency around the time of 
the density jump. 

Our conclusions are very different from those of previous 
works, which predicted a significant optical rebrightening, and 
rather sharp features in the afterglow light curve. The main cause 
for this difference is our careful consideration of both the effect of 
the reverse shock on the dynamics (which we find cannot be ne- 
glected even when a ~ 2), and the arrival time of the photons to the 
observer from different parts of the emitting regions. Both of these 
effects tend to smoothen the light curve significantly. 

Finally, we considered the implications of our results for the 
origin of the fluctuations in the highly variable light curves of four 
GRBs (3 long-soft GRBs and one short-hard GRB). We find that 
density variations are unlikely to be the source of the fluctuations 
in any of these bursts. 
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APPENDIX A: A ONE DIMENSIONAL SPECIAL 
RELATIVISTIC HYDRODYNAMICS AND RADL4TION 
CODE 

In order to obtain accurate light curves while relying on a minimal 
number of approximations we use a one dimensional special rela- 
tivistic hydrodynamic code and combine it with an optically thin 
synchrot ron radiation module. T his is the same code that was used 
before in lNakar&Pirar]<2004 . We use a one dimensional hydro- 
dynamic code that was generously provided to us by Re'em Sari 
and Shiho Kobayashi. It is a Lagrangian code based on a second- 
order Gudanov method with an exa ct ultra-relativistic Riem ann 
solver and it is de scribed and used in lKobavashi et all il999h and 
iKobavashi &"SMiU200QV On top of this code we have constructed 
a module that calculates the resulting optically thin synchrotron ra- 
diation. The code does not include the synchrotron self-absorption 
or synchrotron self-Compton processes. The effect of the radiation 
on the hydrodynamics is neglected. Below we describe the physics 
that is included in the radiation module. 

Having the full hydrodynamic evolution of the fluid (from the 
hydrodynamic code) we first identify the time steps in which a 
given fluid element is shocked by finding episodes of increase in 
its entropy. The same fluid element can be shocked many times. 
Once a fluid element is shocked all its electrons are assumed to be 
instantly accelerated into a power-law energy distribution with an 



index p, dN /d'ye oc 7^'' for 7,, > 7niin. The energy in the electrons is 
taken as a constant fraction, e,., of the internal energy, and this con- 
dition determines 7„,i„ (it is assumed that p > 2). From this point 
on, and until the same fluid element is shocked again, the electron 
energy distribution decouples from the internal energy and evolves 
through radiative cooling and adiabatic cooling or heating {PdV 
work). The magnetic energy in each fluid element is taken to be a 
constant fraction, eg, of the internal energy at all times. 

One of the main difficulties in calculating the synchrotron ra- 
diation at high frequencies is the short cooling time, which may be 
much shorter than the hydrodynamic time steps. In order to over- 
come this difficulty, we calculate the radiation during any hydro- 
dynamic time step analytically, in the following way. Immediately 
after a fluid element crosses a shock, its initial electron energy dis- 
tribution is taken to be a power-law (with index p) between 7„,i„ and 
7max = 00. The total emissivity of the fluid element at a given fre- 
quency in its own rest frame, during a time step, is obtained by inte- 
grating the spectral power of individual electrons over the evolving 
electron distribution, where each electron is tagged by the value of 
its initial Lorentz factor. The emission of each electron is obtained 
by time integration over its instantaneous emissivity', which in turn 
depends on the evolution of its Lorentz factor (and thus on its initial 
Lorentz factor) during the time step. This evolution is calculated by 
considering its radiative losses and its adiabatic cooling or heating. 
In particular, we calculate the evolution of an electron with initial 
Lorentz factor 7mi„ and of an electron with initial Lorentz factor 
7max- Their values at the end of the time step are taken as the initial 
values for the next step in which the initial distribution of electrons 
Lorentz factors is taken again as a power-law between the new val- 
ues of 7min and 7max. From the point where 7max becomes compara- 
ble to 7mi„ (within a factor of 2) the electron energy distribution is 
taken as a delta function. 

Since we use a one dimensional code in spherical coordinates, 
which explicitly assumes spherical symmetry, each fluid element 
represents a thin spherical shell. Once the rest frame spectral power 
of a fluid element is calculated, we integrate over the contribution 
of this shell to the observed flux at a given observer time and ob- 
server frequency. This calculation takes into account the appropri- 
ate Lorentz transformation of the radiation field and photon arrival 
time from each point along the shell. 



' The synchrotron spectral p ower [erg/Hz/sec] of each el ectron is approx- 
imated in the usual manner (Rybicki & Lightmanlll986l) : oc Siusyn — 
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